Resonance Zones and Lobe Volumes for Volume-Preserving Maps 
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We study exact, volume-preserving difFeomorphisms that have heterochnic connections be- 
tween a pair of normally hyperbolic invariant manifolds. We develop a general theory of lobes, 
showing that the lobe volume is given by an integral of a generating form over the primary 
intersection, a subset of the heterochnic orbits. Our definition reproduces the classical action 
formula in the planar, twist map case. For perturbations from a heterochnic connection, the 
lobe volume is shown to reduce, to lowest order, to a suitable integral of a Melnikov function. 

1 Introduction 

The computation of the volume of incoming and exit sets for a "nearly invariant" region is the first 
step in the development of a dynamical theory of transport. For area-preserving maps, it is common 
for these regions to be bounded by segments of the stable and unstable manifolds hyperbolic 
invariant sets, typically periodic orbits or cantori (Aubry-Mather sets) |MMP84l IRKW88] . In this 
case the resulting set is a "resonance zone" |MMP87l IEas91] , and its exit and entrance sets are 
"lobes" of the "turnstile." When the splitting between the manifolds is small, the resonance zone 
is nearly invariant. 

The volume of a lobe is a not only a measure of the separation of the stable and unstable 
manifolds, it is also the flux of trajectories escaping from the resonance zone. This flux provides 
an estimate for the escape time from the resonance. It is known that the average exit time from a 
region is equal to the ratio of accessible volume of the region (the fraction of the volume that can 
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be reached by orbits from the outside) to the flux |Mac941 IMei97j . Though the accessible volume is 
diflttcult to compute it is certainly bounded by the total volume of the region, so the flux provides 
an upper bound on the average exit time. Moreover, if the flux goes to zero, but the accessible 
volume does not, then the average exit time must go to inflnity. 

In this paper, we discuss the construction of lobes for resonance zones of volume-preserving 
maps. Our goal is to generalize the results of |MMP84] that provide formulas for the lobe areas 
for the two-dimensional maps and of |Mac94] for the case of three-dimensional, incompressible 
vector flelds. Indeed, one of the three open problems posed by MacKay at the end of |Mac94j is to 
generalize his flow results to the case of maps. 

The theory of |Mac94j applies to exact-symplectic maps. Recall that a map is symplectic when 
there is a closed two form, w, (for example oj = dq f\ dp) that is preserved by /: f*uj = wj^ A map 
is exact-symplectic when uj is exact, lu = —dv (for example v = pdq) and there exists a function S 
defined by 

f*u-u = dS. (1) 

An often studied case is that of twist maps, which have Lagrangian generating functions — the 
discrete analogues of the Lagrangians for differential equations [Mei92] . 

Remarkably, the generating function ([T| provides a way to compute the area of a lobe (or of the 
resonance zone itself) in terms of the action, the formal sum of S along an orbit jMMP84l IMMP87] . 
The result is that a two-dimensional integral over a lobe is reduced to a zero-dimensional integral, 
the difference between the actions of the orbits homoclinic to the hyperbolic invariant set. Thus to 
compute the lobe area one only needs to find these homoclinic orbits and carry out the sum; this is 
considerably easier than constructing the entire lobe boundary and computing the two-dimensional 
integral. Moreover, since the action is stationary on the orbits, its computation is second-order 
accurate. 

For the case of incompressible vector fields, a similar result also holds. On the manifold M = M^, 
when V-ii = there exists a vector potential A so that u = V x A. Here A is more properly thought 
of as a one- form f3 = A-dl and it gives rise to a variational principle for orbits of the fiow. Moreover, 
MacKay has shown that the one-form (3 can be used to compute the fiux through two-dimensional 
surfaces jMac94] . More generally, when $7 is an exact volume form and u is an exact incompressible 
vector field, then the curl relationship generalizes to the statement that i^O, = d(3 is exact and (3 
provides the generating form (see {7.2 and |LM08] ) as well as a variational principle jGM03j . 

In this paper we study exact volume-preserving maps. For such maps there is also an ana- 
logue of ([T]) in which the generator S becomes a differential form A, see ^ Such forms were 
implicitly described by |Car04j in the context of constructing implicit generating functions for 
volume-preserving maps |LM08] . We will primarily study the three-dimensional case where the 
generator A is a one-form. 



^We recall the notation for the puUback, /*, and similar concepts in the Appendix. 
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We will construct a resonance zone based on a pair of normally hyperbolic invariant manifolds, 
say A and B. For the three-dimensional case, the natural objects are periodic orbits and invariant 
circles. The boundary of the resonance will consist of pieces of the codimension-one stable PF'^(A) 
and unstable W'^{B) manifolds of these invariant sets. The exit and incoming sets are obtained by 
iterating the resonance, see Q 

Unlike the two-dimensional case, the boundary of a resonance does not always consist solely of 
pieces of stable and unstable manifolds. Indeed for this to happen, the manifolds would have to 



intersect on a "proper boundary", see [4.2. However, we have observed that the set of heteroclinic 
intersections commonly contains components that are themselves bi-asymptotic to the invariant 
sets; these preclude the existence of heteroclinic proper boundaries |LMOO| . A similar phenomena 
occurs for three-dimensional volume-preserving flows: when the two-dimensional manifolds of saddle 
equilibria intersect, they typically do so only along a few "primary" heteroclinic orbits |Bro8H 
IHol84| IMac94j . To remedy this, a "cap" must be used to complete the resonance boundary. Even 
though the cap is — to a large extent — arbitrary, this choice does not change the volume of the exit 
and incoming lobes of the resonance. 

Our main result is Thm. |4] in f|5| which states that the volume of the exit and incoming lobes 
is given by 

fcg^A JV(A,B) 

Here the set V{A,B) = [J^^z f^iv) consists of the primary heteroclinic intersections of the man- 
ifolds VF^(74) and W^{B), and r] is the restriction of V{A,B) to a fundamental domain. Thus to 
compute the volume of n-dimensional lobes one needs only to do (n — 2)-dimensional integrals along 
submanifolds of heteroclinic intersections. 

For the case of nearly-integrable systems, a widely used technique for detecting such heteroclinic 
intersections is the Poincare-Melnikov method |Mel63| IHM821 ILMRROS] . Indeed, this method 
detects — at first order — precisely the primary intersections. 

The classical Melnikov function computes the rate at which the distance between the manifolds 
changes with a perturbation, say 6. For two-dimensional maps, an integral of the Melnikov function 
between a pair of zeros gives the rate of change of the flux with 5 |MM88] : a similar result also 
holds for incompressible vector fields jMac94] . We will show in ^that this same results holds for 
volume-preserving maps. 

We conclude by presenting several examples and applications in ^ 
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2 Exact Volume-preserving maps 



2.1 Definition 

Let M be an n- dimensional manifold. A volume form i7 is a nondegenerate n-form on M. A map 
/ : M ^ M on an n-dimensional manifold preserves the volume form 0, if 



For example if M = M", and 0, = dxi A dx2 A ... A dxn, then / is volume preserving when its 
Jacobian has unit determinant everywhere, det(D/) = 1. Suppose now that 0, is exact, i.e., there 
exists an {n — l)-form a such that $7 = da. By analogy with the symplectic case ([T]), we can also 
define exact volume-preserving maps. 

Definition 1 (Exact Volume Preserving). A dijfeomorphism f : M ^ M is exact-volume pre- 
serving if there exists an (n — l)-form a such that da = 0, and a generating (n — 2)-form X such 
that 



It is clear that if / is exact-volume preserving, then / is also. Moreover, if / = ffi o 52 is the 
composition of exact volume-preserving maps with generating forms Ai and A2, respectively, then 
since {gi o g2)* = g^ 05*, 



Thus / is exact- volume preserving with \ = g2\i + \2 ■ Consequently, the set of exact volume- 
preserving maps is a subgroup of the group of volume-preserving diffeomorphisms, we denote it 
DiffQ(M). We discuss some of the consequences of exactness in |LM08j . 

As an application that we will need later, consider the composition of an exact volume-preserving 
map with itself. In this case we can use a telescoping sum to conclude 



f*n = n. 



f*a — a = d\. 



(2) 



fa -a = gligla -a)+gla-a = d{gl\i + A2) . 




i=o j=o j=o 



Thus /" is exact-volume preserving with the form 



n-l 




(3) 



i=o 
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2.2 Geometrical Implications 

Exact volume-preserving maps arise naturally in the context of perturbations from integrable maps. 
For example suppose that / is a map on the phase space M = T*^ x M with volume form = 
dz A dOi A ... A dOd- This form is exact with a = z A dOi A ... A dOd- One integrable map on M is 

f[e,z) = [e + p{z),z) , (4) 

with the rotation vector p : M ^ T''. This map is also exact volume preserving: if we define the 
vector field W{9, z) = (0, / zp'dz), then f*a — a = dX where 

A = iwdOi A ... A dOd . 

The invariant tori of this map are "rotational" tori, and as we argue next, any volume preserving 
map on M with a rotational invariant torus must be exact. 

A rotational torus on T'^ x R is a d-dimensional torus that is homotopic to the zero section 
{{0,0) : 6 G T'^}. The net flux crossing a rotational torus T is the difference between the signed 
volume below f(T) and that below T: 

Flux(T) = j f*a-a. (5) 

Note that if / is volume preserving, and T and T are any two rotational tori then the volume 
contained between them, AV = a — a, is invariant. This implies that Flux(T) is independent 
of the choice of torus. 

Thus, if Flux(T) = then /(T) n T 7^ for any rotational torus; this intersection property is 
used in the generalization of KAM theory for exact volume- preserving maps |Xia92] . Conversely, 
if / has an invariant rotational torus then Flux(T) = by definition. Consequently if / has a 
rotational invariant torus then it is necessarily exact. 

Exactness can be used to simplify the computation of the volume of certain regions. Suppose 
that M is any orientable manifold, and 7^ is a region whose boundary can be decomposed into 
pieces that are related by iteration, S and and such that S H f{S) = C is an invariant 

codimension-two submanifold. We assign orientations to TZ and 5, and by iteration to f{S); thus 

dn = f{s) - s , 

and dS = C. Then, when / is exact, Q immediately gives 

Vol(7^) = / J7 = / f*a-a= [ X. 
Jn Js Jc 

Thus we can compute the volume of TZ simply by integrating over C. This formula is closely related 
to those that we will use to compute the volume of exit and incoming sets in ^ 
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2.3 Examples 

Since every closed form on M = M" is exact, the volume form 0, = dxi A . . . A dxn is exact with, for 
example, a = xidx2 A ... A dxn- If / is volume preserving, then /* a — a is a closed (n — l)-form, 
but since every closed form on M" is exact, / is exact. 

However, exactness is not automatic on more general manifolds. For example, a perturbation 
of Q is the one- action map / on M = T'^, xM 

f[e,z) = [e + p[z + Fm.z + Fm , (6) 

for a "force" F : T*^ ^ M is always volume-preserving, but is only exact when the form FdOiA. . .AdOa 
is exact |LM08] . or equivalently when 

/ Fd9i A...Aded = 0. 

For example, if d = 2, p{z) = (z, z^) and F = a cos 6i + b cos 62 + cos(0i -|- 62) then ^ is generated 

by 

A = Qz^ + asin^i^ (i^2 - + bsm92 + csin(^i + ^2)^ dOi , 

where Z = z + F{9). 

In |LM08j we showed that any exact-symplectic map ([T| is also exact-volume preserving. For 
example, reinterpreting ([6| as a map on {9, z) G x M^, then / becomes a "generalized standard" 
or "Froeschle" map |Fro72j . It is symplectic with the two-form lo = Y^i^-Gi ^ dzi provided that 
p{z) = VK{z) for "kinetic energy" K and F = — W{9) for "potential energy" V . Letting v = z-d9, 
then / is also exact-symplectic and has the generator 

S{9, z) = Z- VK{Z) - K{Z) - V{9) , 



which is equivalent to the implicit, Lagrangian generating function [Mei92j . For example when 

1 1 



p{z) = z, then K{z) = 4|zp and the generator becomes 



S{9,z) = ]^\Z\''-V{9). 

Now setting a = A w, then da = O = d9i A d92 A dzi A dz2 and we find that the generator (|2 
for the generalized standard map is the two-form 

X=^Si9,z)uj. 
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3 Normally hyperbolic invariant manifolds 



Resonance zones are most naturally associated with normally hyperbolic invariant manifolds; in the 
simplest case, with a hyperbolic periodic orbit. In addition to periodic orbits, we will also consider 
normally hyperbolic invariant circles, but many of our results apply more generally. Here we recall 
the definition of normal hyperbolicity, and prove a lemma about convergence of forms on the stable 
manifolds of such sets. 

Suppose and / : M — > M is a diffeomorphism on an n-dimensional smooth manifold M and A 
is a compact invariant set of /. We recall one standard definition of normal hyperbolicity. 

Definition 2 (Normal Hyperbolicity jPdlLSOSj ). A compact invariant set A is r- normally hyper- 
bolic for r £ N, if there exists an invariant splitting of the tangent bundle TaM = (B © TA, 
a Riemannian structure, and positive constants C, A and n, such that for all a £ A, 

a) < A < /i"'' < 1; 

b) \\Df\a)v\\ < CA"||w||, for all v £ El and t £ N; 

c) \\Df-\a)v\\ < CA"||u||, for all v G E^^ and t G N; 

d) \\Df\a)v\\ < Cn^''^\\v\\, for all v £ TaA and t E Z. 

The stable manifold theorem applies to r-normally hyperbolic sets: there exist C immersed 
submanifolds WiAJ) = WiA), tangent to E"" ®TA at A, and VF"(yl,/) = VF"(yl), tangent to 
E'^eTA at A |HPS77aj . 

We next prove a lemma that we will need in the following sections to show that some of the 
sums converge. This lemma applies to differential £-forms on a compact part of a stable or unstable 
manifold of a normally hyperbolic invariant set. 

Let A^(M) denote the linear space of £-forms on M. For any compact set P C M, there is a 
natural norm on A^(P) given by: 

\\uj\\p = sup{\uJp{vi,V2, . . ■,Vi)\ -.pGP, Vi £ TpP, \\vi\\ = 1, i = l,...,i} . 

The main point is that for any p £ P and any set of vectors Vi £ TpP, 

\uJx{vi,V2, . . .,Vi)\ < \\uj\\p \\vi\\ \\V2\\ ■ ■ ■ \\Vi\\ . 

Definition 3 (Regular Form). A differential l-form uo defined on W^{A) is regular if 

lim {f*fLO = . 

fc— >oo 

Similarly if uo is defined on W^^{A) then it is regular if 

lim {f*fuj = . 

fc^— oo 
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In ^ and ^ we will use the following simple consequences of this definition. 
Lemma 1. If lo is regular on W^{A), then 

a) if P d W^{A) is a compact suhmanifold of dimension i then 

lim / = ; 

t^OO Jp 

b) IfY is a vector field on W^{A) and Q C VF'^(^) is submanifold of dimension £ — 1 then 

lim / iy(/*)*a; = 0. 

An £-form is automatically regular if £ is large enough so that every independent set of i vectors 
Vi G TW^{A) is guaranteed to contain sufficiently many vectors in the stable, as opposed to center, 
directions. 

Lemma 2. If A is an r-normally hyperbolic invariant set with dim(^) = ua, then every differential 
£-form with £ > and 

£>nA(l + -] (7) 



is regular on any compact subset P C Ty^(A). 

Proof. By hypothesis there is a splitting TaM = (B © TA and positive constants C, A and 
/i as in Defn. [2j The stable manifold theorem implies that there exists a neighborhood AA of A in 
W^{A) and coordinates cf) = {a, s) : J\f ^ A x M"", where = dim(£''^), such that s(^) = if and 
only if ^ G j4, and such that / = o / o takes the form 

/(a, s) = {g{a) , L{a)s + r{a, s)) 

with r(a, 0) = 0, d2r{a, 0) = and L{a) is a matrix. In addition, there is an induced Riemannian 
structure on M such that Df satisfies the conditions of normal hyperbolicity on the zero section, 
for some constants C, A and fi, such that A<A<1, 1 < < and fTX < 1. 

Thus if P is a compact submanifold of VF'^(A), the normal hyperbolicity conditions imply that 
for all t G N there exists mi , m2 G N such that 

||(/*)*a;||p<C'(/i™^A™^)*||u;||p. 

where mi + m2 = £ and mi < n^. The r*^ power of the contraction factor in this equation satisfies 
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and ([7| implies that 

— mi = {i — mi)r — mi > {i — n^)r — ua > . 

Thus both integers mi and m2r — mi are nonnegative. If at least one of these integers is positive, 
then /i^iA™^ < 1. If mi > then all is well. Alternatively, if mi = 0, then m2r — mi = ir, which 
is positive since £ > hy hypothesis. Consequently, in either case, ||(/")*u;||p — > as required. 

□ 

Remark. It is important to notice that we have to include the manifold A in otherwise we 

could have invariant forms on iy^(j4) that don't satisfy the conclusion of Lem. [2] For instance, 
in |LM03j . an invariant n — 1 form is constructed on \ A when ^ is a normally hyperbolic 

invariant circle. We also want to point out that similar computations appear in [DdlLSOS] and 
(CMLQI . 

If UA = 0, then Lem. [2] implies that all ^- forms with £ > 1 are regular, but if = 1, then £ >2. 
We will often apply the lemma to the case £ = n — 1, in which case an £-form always is regular on 
a one-normally hyperbolic invariant manifold when riA < (n — l)/2. 

4 Measuring Transport 

In dynamical systems, transport is the study of the motion of collections of trajectories from one 
region of phase space to another. A natural choice for a region is a subset IZ <Z M that is almost 
invariant in the sense that it consists of a set of points whose orbits belong to IZ for a long time. One 
way to construct such a nearly invariant region is to form its boundary dlZ (as much as possible) 
from the invariant manifolds of a pair of normally hyperbolic invariant manifolds; these boundaries 
are called partial harriers^ 

For the two-dimensional case this leads to the construction of resonance zones bounded by 
the invariant manifolds of a periodic orbit or cantorus |MMP84l IEas91] . For a saddle fixed point, 
o, it is natural to base the selection of the exit and incoming sets on primary homoclinic points 
p,q ^ W^{a)r]W^{a), see Fig.[TJ In this case the boundary of the resonance is a single partial barrier 
formed from segments of the stable and unstable manifolds from a to f{p). In other cases, like 
maps on the cylinder, the resonance zones associated with rotational periodic orbits have upper and 
lower boundaries, each formed from different branches of stable and unstable manifolds |MMP87] . 

In this paper we will generalize this construction to the volume-preserving case. To do this, we 
will not be able to assume that the lobes are bounded entirely by subsets of the invariant manifolds. 

In order to restrict the topological possibilities to something manageable, we will usually assume 
that the phase space M is at most three-dimensional and consider only the simplest hyperbolic 
invariant sets: fixed points and invariant circles. 
^ Another is to use an isolating block [Eas91| 
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Figure 1: Standard construction of a two-dimensional resonance zone TZ for a saddle point a, based 
on a primary homoclinic point f(p) G W^{a) D W^{a). The boundary of the resonance zone TZ is 
W^, J a) U Wf, -.(a) and TZ has exit set £ and incoming set Z. 



4.1 Incoming and Exit Sets 

Any region TZ C M has incoming and exit sets. The exit set £ is the set of points in TZ that leave 
in one step: they have forward exit time one. Similarly, the incoming set X consists of all points 
not in TZ whose images land in 7^0 

Definition 4 (Incoming and exit sets). The exit and incoming sets for a region TZ are 

£ = TZ\f-HTZ), 
I=r\TZ)\TZ, 

respectively. 

Suppose now that / preserves a measure fi, and that fi{TZ) < oo. Since fi{A \ B) = fJ,{A) — 
fj.{Ar\B), then 

fi{£) = fi{TZ)-fi{TZnf-Hn)) , 

MX) = Mr 'W) - Mr ' W n 7^) . 



^ In [Mei97] . the entry set was defined to be the image of the incoming set and is a subset of TZ. Here it seems 
more convenient to use the preimage of the entry set. 
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My measure preservation, ij,{TZ) = n{f^^{TZ)), so that the volume of the exit and incoming sets are 
equal: /x(X) = n{£). 

The incoming and exit sets for a region are often called "lobes" |RKW88] , and their union is a 
"turnstile" |MMP84] . The volume of a lobe — either exit or incoming set — is a simple measure of 
how fast the points move from TZ to its complement. In other words, the lobe volume is a measure 
of the degree to which a set is invariant. Roughly speaking, if one selects a point at random in TZ, 
then it has probability p = ^{8)/ ^ilZ) of landing in £■ at each iteration, and so its expected escape 
time is 1 /p. This can be made precise |Mei97j : the average time for a trajectory that begins in X 
to escape is 

where TZacc is the accessible portion of TZ: the part that can be reached by trajectories that begin 
outside TZ. To obtain more details of the distribution of exit times requires an understanding of 
the decomposition of Z into regions of fixed exit time |RKW88| IMei97] . 

4.2 Fundamental Domains 

To construct the boundaries of a resonance zone, we will use the concepts of fundamental domain 
and primary intersection of the stable and unstable manifolds of a hyperbolic invariant set |LMOO] . 

In order to form a resonance zone based upon a pair of normally hyperbolic invariant manifolds 
A and B, they must have co dimension-one invariant stable and unstable manifolds. To be concrete, 
we will suppose that 

codim(VF"(.4)) = codim(VF'^(B)) = 1 . (9) 

If cr C H^''(^) is the boundary of a subset of M^'^(A) that contains A, then we denote the unstable 
manifold starting at a by W^{A): it is the closed subset of the local stable manifold of A bounded 
by cr: 

a = dW^M) ■ 

Similarly, W^{B) is the unstable manifold up to 7 when 7 = dW^{B)\ however, in this case it is 
convenient to assume that this submanifold is open. Though this definition is not not symmetric, 
the asymmetry is useful to simplify some proofs. 

For our purposes, it will be important that the boundaries of these local manifolds are chosen 
to be proper. Recall that a neighborhood of an invariant set A is isolating or a trapping region 
\i Ac int(iV) and 

/(cl(iV)) C int(iV) . 

Definition 5 (Proper Boundary). Suppose A and B are compact, normally hyperbolic invariant 
manifolds. A set a C W^[A) is a proper boundary if W^{A) is an isolating neighborhood of A in 
W^[A). Similarly, 7 C W^{B) is proper ifW^[B) is an isolating neighborhood of B for /"^. 
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It is not hard to see that the stable manifold theorem, e.g. |HPS77bj . implies that proper boundaries 
always exist. These local manifolds behave naturally under iteration: 

fiW'M)) = W}^.M) , and f{W^{B)) = Wf^^){B) . 

Given a proper boundary, an invariant manifold can be partitioned into nonoverlapping fundamental 
domains that are related by iteration. 

Definition 6 (Fundamental Domain). The set of fundamental domains ofW^{A) is the 

collection of sets of the form 

where a is any proper boundary. Similarly, the set J-^{B) of fundamental domains in W^{B) is 
the collection of sets of the form 

u^{B)^wf^^^{B)\w:;{B), 

where 7 is any proper boundary. 

Remark. The closure assumptions imply that a C and j C U^, but their images are not. In 
addition, one has dS^ = cr U f{a) and dU-y = 7 U 7(7). 

As an example, consider the two-dimensional map sketched in Fig. [T} In this case, the resonance 
zone TZ is bounded by U ^/(p)!*^) the fundamental domains Sp{a) and Up{a) form the 

boundary of the region £ U I, the turnstile for TZ. For a three-dimensional map, a fundamental 
domain is an annulus bounded by a proper boundary 7 and its image 7(7). 

A consequence of the definition is that the image of a fundamental domain is also a fundamental 
domain and that f^{Sa) = •Sfk(^„^ and f^iU-y) = ^fk(-y) for any k £ Z. Moreover, stable and unstable 
manifolds can be decomposed as the disjoint union of fundamental domains: 

(10) 

w:;{B)\B = [jup(^y){B). 

t<o 

Consequently, the topology of the intersections of stable and unstable manifolds can be studied 
by restricting to appropriate fundamental domains |LMOO| ILM03] . We will use the fundamental 
domains to construct the incoming and exit sets for a resonance zone. 

4.3 Primary Intersections 

As above, we continue to assume that / has a pair of normally hyperbolic invariant sets A and B. 
In addition, the stable and unstable manifolds ^^^(A) and W^{B) are co dimension-one, orientable 
manifolds as in (|9]). In addition, we now assume that there is a heteroclinic intersection 

W''{A)nW''{B)^$. (11) 
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The set of heteroclinic intersections is typically staggeringly complex; in this section we pick out 
the first or primary intersection to use in the construction of a partial barrier from these manifolds. 

For example, consider a pair of saddle fixed points a and b, of a two-dimensional map. A point 
rj £ W^{a) n W^{b) is a primary intersection point (p.i.p.) if (recall W^{b) is open and W^{a) is 
closed) 

T^,^(a)nW,"(6) = 0; (12) 
that is, the manifolds up to r] intersect only at their boundary, |Wig92 . For example, the points p 



and q are primary homoclinic points in Fig. [T] but the point r is not since the set W^(a) Pi Wll^{b) 
contains five points. This definition will not work in higher dimensions, since the components of 
the set 14/^'^(^) n W^{B) need not be proper boundaries even when the intersections are transverse. 
Instead, we generalize as follows: 

Definition 7 (Primary Intersection). A point rj £ 1^*^(71) riW^{B) is a primary intersection point 
if there exist proper boundaries 7 and a such that 

V G w'^iA) n wf^^^iB) , but w'^iA) n w:;{b) = . (13) 

The set of primary intersections is denoted V{A,B); it is an invariant set. 



For one-dimensional manifolds, this gives the same set as (12). Primary intersections can also be 
defined in terms of fundamental domains. Indeed, since = SaD W^^^-^ and = U W}^ 



by (10), then if is a primary intersection 

V G {s,{A) u W}^^^{A)) n {u,{B) u w:;{B)) = s^{A)nu^{B) , 



since, by (13) S„{A) n W^(S) = U^{B) n VF'|(^)(^) = 0. In addition, (10) and (13) imply that all 
forward images of the stable fundamental domain are disjoint from the unstable one: 5jt(o-)(^) ^ 
U^{B) = for t > 0. This means that the intersection index, 

k{SM) = sup{t G Z : f\S) n / 0} , (14) 

is zero. Consequently, an alternative characterization of the primary intersection set is |LMOO| 

V{A,B) = {Sr\U: k{S,U) = 0,5 G r^{A)M G J^"(S)} . 

We will henceforth make the assumption that the co dimension-one manifolds VK*^(A) and W^{B) 
have transversal primary intersections (the manifolds may still have tangencies elsewhere). In this 
case, the components of V{A^ B) are codimension-two submanifolds that never cross. For example, 
in the three-dimensional case, these components must be closed loops or else are curves that are 
asymptotic to the invariant sets A and B. 



13 



If the heteroclinic intersections appear as the result of spUtting of separatrices, then typically it 
is possible to apply Melnikov's method. Recall that there is a correspondence between the zeroes of 
the Melnikov function and heteroclinic intersections. If all the zeroes are simple then they continue 
precisely to the set of primary intersections. 

Since each fundamental domain generates the entire manifold, we can restrict attention to the 
set of primary intersections in a particular fundamental domain, say VCiU^. A fundamental domain 
on a two-dimensional manifold is an annulus, but under the natural identification 7 ~ /(t) of its 
boundaries, it can be thought of as a torus U^, see Fig. [2] 

If ?? is a component of V then, under this identification, its orbit /*(??) is equivalent to a closed 
loop. Indeed, if t] is contained in the interior of cl(U^), then it must be a loop since the intersections 
are assumed transverse. On the other hand, if there is an intersection point p = r] Ci'j, then since 
V is invariant, the curve /(//) G V and intersects f^j) at the point /(p). Thus the components r] 
and f{rj) are joined by the natural identification. 

Continuing this implies that the full orbit of p lies on 7 under the identification. This orbit 
must be finite, since it would otherwise have limit points, violating transversality. Thus, the orbit 
of r] becomes a closed loop fj on U^. 




Figure 2: Identifying the boundaries of a fundamental domain gives a torus U and primary inter- 
section loop fj. 

Therefore, on the torus U-y, a primary intersection loop can be labeled by its homology class, 
{m,n) G 1? . A loop that is homotopic to 7 will be said to have class (0,1), while loops that 
correspond to intersections ry that are asymptotic to A and B will have class (m, n) with m 7^ 0. 
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We previously used this classification to discuss bifurcations of V that occur when the manifolds 
develop tangencies as a parameter is varied |LMOO| ILM03| ILRROSj . 



4.4 Resonance Zones 

In this section we will establish the basic assumptions to construct resonance zones in terms of 
codimension one stable and unstable manifolds of normally hyperbolic invariant manifolds. 

The geometry of resonance zones based on arbitrary normally hyperbolic invariant sets could be 
quite complicated. In order to gain some intuition, we start by describing the case that / is a map on 
a three-dimensional manifold M and assume that A and B are hyperbolic fixed points or invariant 
circles. Motivated by this discussion, we will then propose a set of geometrical assumptions to 
define a partial barrier for a resonance zone. 

Three typical resonance zones IZ for this case are sketched in Fig. |3j For example, if A and B 
are fixed points, IZ \s a. ball that is, roughly speaking, bounded by W^{A) and W^{B) for some 
proper boundaries a and 7. The sketch corresponds to the "integrable" case when these manifolds 
coincide forming a saddle connection from Ato B. More generally the boundary of IZ will be made 
from a partial barrier V that is constructed from W^{A) and W^{B) plus a "cap" (see below). 




Figure 3: Resonance zones in the integrable approximation for a pair of fixed points, a single 
invariant circle or a pair of invariant circles. 

When A = B \s a hyperbolic invariant circle and Ty^(yl) H 7^ 0, the resonance zone is 

a solid torus. This is also true when A ^ B, however, the boundary in this case is obtained from 
a pair of partial barriers, P constructed from VF'^(A) and W^{B) and V constructed from 
and W^{B) (which we also assume intersect); as before the simple case where the partial barriers 



15 



are saddle connections is sketched in Fig. [3} If j4 is a fixed point and B is an invariant circle, then 
TZ will typically be a ball, but a second fixed point or invariant circle will be needed to complete 
the resonance, and its boundary will consist of at least three partial barriers. 
In the general case we require the following geometrical hypotheses. 

(HI) A and B are normally hyperbolic invariant manifolds of dimension at most ^^y^.; 

(H2) l^'^(A) and W^{B) are orientable, codimension-one submanifolds; and 

(H3) the set of primary intersections V{A, B) C VF'^(^) n W^{B) is transverse. 

To construct partial barrier P, begin by selecting a pair of proper boundaries a and 7 such that 
the associated fundamental domains S = Sa{A) and U = U-f{B) are in "standard position", i.e., 
such that 

(H4) k{S,U) = and the set i] = S ^^U C V a. neal|^ submanifold of S and U: 

dTjCdSndU . (15) 

This condition is sketched in Fig. |4j 

The simplest case corresponds to the primary intersection set containing a proper boundary, 
for then we can choose 7 = a C P. In this case (H4) is automatically satisfied since d{S CiU) = 0. 
Recall that this is what is typically done for the two-dimensional case, as shown in Fig. [l] However 
V does not always contain a proper boundary, i.e., a curve with homology (0, 1); in particular, we 
commonly observe that V contains families of curves that spiral asymptotically from B to A, as 
sketched in Fig. [2] |LM001 ILM03j . This also must occur when the map is a Poincare map of an 
autonomous flow, since heteroclinic points of the map lie on heteroclinic orbits of the flow |Mac94j . 

When resonance zone cannot have a boundary that consists solely of pieces of stable and unstable 
manifolds we must add a cap C to construct a partial barrier, 

V = W^{A) U W^{B) U C , (16) 

that will be one of the boundaries of TZ, see Fig. |4j We assume that is it possible to choose the cap 
C so that 

(H5) C is a codimension-one submanifold with boundary 

9C = 7 U a ; (17) 

(H6) C n W^{A) = cj and C n d{W^{B)) = 7. 

A submanifold with boundary V C W is neat if it is closed in W and its boundary dV is contained in the 
boundary dW of W. 
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Figure 4: Partial barrier constructed from W^{A), W^{B), and a cap C. 



For the three-dimensional case the partial barrier under these assumptions is topologically a 
sphere when A and B are fixed points, like that sketched in Fig. [3| When A = B is an invariant 
circle then P is a torus, and when A ^ B are invariant circles, D is an annulus bounded by the 
circles. 

In the exceptional case that the set of primary intersections includes a proper boundary, we 



choose a = J C V. In this case, C = 0, and the partial barrier is still given by (16). 



5 Lobe Volume 

If the map is exact area-preserving, the lobe volume depends only upon the orbit of the manifolds, 
and can be computed using the generator S" of ([T]) |MMP84| IMMP871 lEasQlj . In this section we 
will show that the computation of lobe volume for the exact volume-preserving case will reduce to 
the integral of the one-form A of Q along the primary intersection curves. 



Before proceeding, in ^5.2, to obtain the formula for lobe volumes, we first derive an iterative 



formula relating the surface integral of a over a submanifold to those over its images. 
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5.1 Iterative Formula 



Here we obtain a fundamental iteration formula for the integral of a over a codimension-one sub- 
manifold Q that is arbitrary, except that its boundary can be written as dQ = 7 U 7(7) U r/, i.e., as 
the union of two types of pieces: a manifold 7 and its image, and the remaining part t] that is not 
related by iteration. This situation is sketched in Fig. [5| 




Figure 5: ^ is a submanifold such that a piece 7 of its boundary is related through / with another 
piece with opposite orientation, — /~^(7). In the diagram, an orientation of Q induces an orientation 
on the boundary, and so 9^ = 7 — f^^i'j) + rj. 



Lemma 3. Suppose f : M ^ M is an exact volume-preserving diffeomorphism, and Q is 
an oriented codimension-one submanifold with boundary dG = j — f{^) + rj, where 7 and rj are 
codimension-two submanifolds. Then, for any t £ N 




Proof. The composition formula ([3]) implies that 

t-i 

{f'ra-a = J2ifrdX. 

k=0 

Integration of this relation over the region G and applying Stokes' theorem gives 

/ (/*r«- I ^ = Y.( (/')*A. 

The integrals over the images of the boundary curves 7 are a telescoping sum, so that 
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Upon noting that, for example, 

'if 



a = a 
5 JfHG) 



we see that the equation above is just the first hne of (18) rearranged. 

The remaining result can be obtained by a similar iteration, but backwards. Note that ^ 
implies a — {f~^)*a = d{{f^^)*X); this can be iterated to give 

-t 

k=-l 



Integrating this relation over Q, as before, and rearranging gives the final line of (18). □ 



5.2 Lobe Volume Formula 



We continue to assume that / obeys (H1)-(H3) and has a partial barrier V, (16), constructed 
from fundamental domains S and U and a cap C that obey (H4)-(H6). The turnstile for D is the 
union of the exit and incoming sets associated with the barrier; it is bounded by the fundamental 
domains, S and U, and the cap C and its image /(C) 

d{£ul) = UUSUCU f{C) . 

Though C is somewhat arbitrary, the lobe volume will be independent of this choice because its 
boundary contains both C and its image. 

Since M^^(A) and W^{B) have codimension one, in general they separate the manifold M. A 
consistent orientation for Ty^(j4) and W^{B) will define an "outside" and an "inside" of the barrier 
D. This in turn will induce an orientation of a and, by iteration, an orientation of f{cr). In this 
case we can write dS = a — f{cr) and dU = 7 — /(t). 

The exit set is the portion of the turnstile where the unstable manifold is outside the stable 
manifold, we will use a + sign to denote this subset. The dividing set between exit and incoming 
lobes is the primary intersection rj = H S~ = r\U~ C V{A, B), which we have assumed is a 



submanifold obeying (15). Taking into account the orientation we write 

S = S+ + S- and U = U+ +U- , 
see Fig. [6} Thus the exit lobe has boundary 

d£ = U+ -S+ +C+ - f{C+) . (19) 
The primary intersection r] also divides the loops a and 7 into pieces that can be labeled it 

cr = cj"*' + and 7 = 7^ + 7" , 
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as sketched in Fig. |6j Consequently, the fundamental domains and the cap have boundaries 

dS^ = ±r? + a± - /(cT±) , 

aZ^± = ±r? + 7^-/(7^) , (20) 



We will compute the volume of the exit lobe (19), but as we discussed in ^4.1 the volume of 
the incoming lobe is the same. 





-f(C) 



Figure 6: Three-dimensional lobes. On the left is illustrated the boundary of the exit and incoming 
lobes formed from a pair of fundamental domains that are subsets of stable and unstable manifolds. 
On the right a pair of caps are added to complete the exit lobe. 



Theorem 4 (Lobe Volume). Suppose that D is a partial barrier for an exact volume-preserving 
map f, obeying (H1)-(H6) and such that the one-form A is regular on W^{A) and W^{B). Then 



the volume of the exit lobe £ with the boundary (19) is 



Yo\{£) = Y,! (/')*A . 



(21) 



Proof. If and are pieces of fundamental domains with boundaries obeying (20 ), then Stokes' 
theorem with (p|, (|19|), and dC^ = — 7^ gives 



U = I a 

£ Jae 



a — I a + a — I a 

u+ Js+ Jc+ Jf{C+) 

a — a — A 

U+ Js+ JdC+ 



(22) 



/ «+ / A— / a — 

Ju+ Js+ i(T^ 



A . 
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By (20 ), Lem. [s] applies to the terms in the last line of (22 1, so that for all t E N, 

r a+ [ \=f^ fu''r>^+ [ «+ / A, 

JU+ Jy+ k^iU Jf''{U+) Jf-Kl+) 

A = -E /(/'rA+ / 



(23) 



Here we have selected the direction of iteration to take advantage of the contraction of the manifolds. 
Since by (HI), dim(A) and dim(i?) are at most Lem. [2] implies that the (n — l)-form a is 
regular on W^A) and W^^B). Thus, since cl(Z^+) C W^^B) and cl(5+) C W^A) are compact 
Lem. [1] implies 



lim 

4— >oo 



a = , lim 



a = . 



Finally, since A is assumed to be regular on M/^'^(A) and W^{B), then 



lim 

t— >oo 



A = 



/-'(7+) 



lim 

t— >oo 



A = 



/*(<^+) 



(24) 



After taking the limit in (23) and substituting back into (22), we find 



P —oo „ oo ^ 



□ 



which is equivalent to (21). This concludes the proof. 
Remark (1). The primary intersection between A and B is given by 

V{A,B)=[jf'{rj). 
Therefore, under the assumptions of Thm. |4] (21 ) can be written 

JS ^Jr? JV{A,B) 

Remark (2). By Lem. [2| the (n — 2)-form A is regular on the (n — l)-dimensional manifolds H^'^(A) 
and W^{B) whenever dim(A), dim(S) < ^ — 1. This includes the case that A and B are fixed 
points and n > 3; however, A is not necessarily regular for n = 3 when A oi B are invariant circles. 
When A = B, the hypothesis that A is regular is probably not needed. It is only the difference 



Jf-t{j+) Jf 



A 



that must limit to zero in (23) to obtain (21). This difference converges exponentially to zero if 
/^*(7^) A and f*{a^) A as t ^ oo so that the domains of integration become identical. This 
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is what happens for n = 2 where A is the zero-form S, ([T|, |MMP84l IMMP87] . and for the examples 
m ^ It would be nice to show that it is always possible to select proper boundaries so that this 
is the case. When A ^ B, it is not obvious whether the requirement that A be regular can always 
be satisfied. Since any closed form can be added to A without changing it may be possible to 
use this freedom to make A regular. One way to do this would be to have X\a = Mb = 0. 

Thm. |4] implies that the volume of an n-dimensional exit set can be computed by integrating 
along the (n — 2)-dimensional primary intersection set. In general, the volume of an n-dimensional 
lobe can be computed by integrating the (n — 2)— form A on the primary intersection V{A,B) 
that is an immersed submanifold of dimension n — 2. Moreover, the flux across a partial barrier 
determined by the manifolds VK'^(A) and W"{B) is independent of the selection of the fundamental 
domains and of the cap C. 

Finally, note that Thm. [4] also applies to the case that the cap C is the empty set. This would 
be the case if V includes a proper loop, for then we can select o" = 7 G P and V = W^{A) U W^{B). 



6 Melnikov Flux 



In this section, we show how our formula (21) for the lobe volume limits to well-known Melnikov 
results for maps with a near saddle-connection. This is well-known for the two-dimensional case 
|MM88j and is implicit in the theory developed in |LMOO] for the three-dimensional case. To prepare 
for the result, we recall some notation and results of |LM03| ILMRR08] . 

Suppose : M — > M is a smooth family of exact volume-preserving diffeomorphisms ^ 



satisfying (H1)-(H3) of ^4.4 for each < |(5| < 60. In addition, suppose that the map /o has a 
saddle connection 

S C W%Ao) n PF"(Bo) . 

between the normally hyperbolic invariant manifolds Aq and Bq like those sketched in Fig. [Sj As 
usual, we assume that dim(E) = n — 1 and dim(Ao)) dim(i3o) < 

By (H3), when 5 > the saddle connection splits into the stable and unstable manifolds 
WiAs) = W^As, fs) and VV^iBs) = W^iBs, fs) of the perturbed invariant sets. The splitting of 
these manifolds is computed, to lowest order, by the classical Melnikov function. 

Specifically, suppose there is a fundamental domain C T, and two diffeomorphisms ■0^'" : C 
E — > M (adapted deformations), such that 

Us = and Ss = V4(-^) 

are fundamental domains of the perturbed manifolds ^^^(jd^) and W^{Bs). When (5 = 0, we can 
take V'o'*^ to be the trivial inclusion -i/'^'" : C S — > M, '0o'^(O = so that 

Un = Sn = T . 
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We notice that the intersection index k(5o,Z^o) = k{J-',J-) = 0, since is a fundamental domain, 
but this need not be true for Ss and Us, unless the deformations were chosen carefully. Nevertheless, 
as in (H3), we assume that the set of primary intersection rjs is transverse for small 5. As we will 
see, the limiting set tiqCi J- is the zero set of a Melnikov function. 

To first order in 6 the only relevant quantities associated with and the deformations tpg'^ are 
their vector fields with respect to 6. The perturbation of the map /o away from (5 = is measured 
by its perturbation vector field 



d_ 
d6 



5=0 



(25) 



for each x G M. The perturbations of the deformations " are similarly defined by the vector 
fields 



d_ 
85 



5=0 

for each £ J', using the assumption that = id on J^. 

Up to vectors tangent to the saddle connection, these vector fields are related by the iterative 
formulas |LMKK08] 



(/o)*z"-^-z"'^ + XoGr.F, 

To see this, note that since S is invariant under /o, the functions ip^ 

adapted deformations, mapping S to W|'". In either case, the function C's = ipj^ o ^5 defines a 
curve in !F, so that the vector -^Cs is in TJ^. Computing this derivative implies that 



(26) 
(27) 

fs ° -05''' o /o~^ are also 

-1 



d_ 
dS 



5=0 



However, 



d_ 
d5 



5=0 



Xo{i) + DUf,\S))Z{f^'{i)) 



Recalling the definition of the pullback and pushfoward (see (61 )) and combining these two results 



gives the pair of formulas (26) and (27). 



Since the maps ip^'^ describe the deformations of the stable and unstable manifolds, the "veloc- 
ity" of the splitting with 6 is given by the deformation vector field 



A(e) = z"(e) - z\o . 



(28) 



The stable and unstable manifolds are thus "flowing" with respect to 6 according to the vector 
fields Z^'^ and they split at a rate A. This flow induces a "flux" with respect to the volume form 
$1 defined by 

^ = 1/^9. . (29) 
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The form $ is the flux form of the deformation A. 

Though the expression for the sphtting velocity A depends on the choice of adapted defor- 
mations, the form is well-defined on T, and is independent of the choice of ■0*'" |LM03] . In 
addition, from the definition (28) with (26) we have 

/o*A - A G ; (30) 

which implies that 

/o (^A^^) = if* A (/o^) = if*A^ = on S , 

i.e., the flux is invariant under /q. 

The flux is our primary measure of the splitting; indeed as we will show below, the integral of ^ 
is related to the lobe volume. Since fs is exact, the net flux across a fundamental domain vanishes. 

Proposition 5 ( [LM03) ). If is a fundamental domain in S then 



$ = . 

T 

The classical Melnikov function is defined relative to a choice of an adapted normal vector field 
n : S ^ TM |LM00j . Recall that vector field is adapted to /o if, for all vector fields Y and all 
points ^ G S, 

/o*(n,y) = (n,/o*y). (31) 

Given an inner product ( , ) and such a vector field, the Melnikov function : S — > M is simply 
the normal component of the deformation, 

Mn = (n,A). (32) 

Since n is a vector field in the algebraic normal bundle of S, is a measure of the speed of 
splitting in the normal direction. Conditions (30) and (31)together imply that is invariant 
under /o, that is, Mn o /o = M^- One common choice for n is the gradient of an invariant of /o (if 
one exists) jLMOO| . but any adapted normal vector field gives a Melnikov function with the same 
zero set. 

The main point of Melnikov theory is that if a point ^ G S is a nondegenerate zero of Mn, the 
stable and unstable manifolds W^{Bs) and ^^^(jd^) intersect transversally near ^ when 6 is small 
enough. The flux form is simply related to the classical Melnikov formula by 

$ = Mnt^n , (33) 

where lo^ is the natural (n — l)-volume form on S induced by n |LM03| ILMRR0"8] . 

LOn = I ^ y in^ ■ (34) 

(n,n) 
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When n is an adapted normal, ujn is invariant under /q: /oWn = on S. 



The iterative formula (|26j) implies jLMn3l ILMRROS] 



A 



oo 

E 

fc=— oo 



(35) 



Indeed, the function A can be thought of as a section on a normal bundle. It was proven in 
|LMRR08] that if one considers both A and the series YlT=-oo ifo)^-^ sections, then they are 
the same function. 

That is, the deformation vector field is, up to tangent vectors, given by the infinite sum above. 
This infinite sum is geometrically convergent since the invariant sets Aq and Bq are normally 



hyperbolic. By (32) and (35) the Melnikov function can be written 

(n,X)o 



oo 

E 

fe=— oo 



fk 

Jo 



(36) 



since, by definition, (n, Y) = for any Y £ TT. 

There is considerable freedom in the selection of the adapted deformations V'^' we may choose 
any functions from T to W^'^ that reduce to the identity when 5 = 0. Given an given adapted 
normal n, there is a natural choice: for each i £ J- choose V'K?) (0) to belong to the intersection 
of the stable (unstable) manifold with the line generated by n(^). Under this choice A = (/in for 
some function : S ^ M to be determined. 

Since A is also a deformation vector field, (35 1 implies that A — A G TT, so that = (n, A), 
as well. From this, we conclude that = {n,n)(j) and therefore A has to be of the form 



A 



(n,n) 



n. 



(37) 



We are now prepared to state relate these old results to Thm. [4] showing that an integral of the 
flux form <I> gives the rate of growth of the lobe volume with 5. 

Theorem 6. Let fs be a family of exact volume-preserving diffeomorphisms with partial barriers 
Ds <ind exit sets £s obeying (H1)-(H6), such that /o has a saddle connection 



S C W\Aq) n VF"(Bo) . 



Lett Mn be a Melnikov function (32) for fs and = T r\ ^ 



be the set of points in the 



fundamental domain T for which the Melnikov function is positive. Then 



d 
dS 



S=0 J£s 



1 



(38) 
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Proof. By assumption fs is exact-volume preserving so there exists a family of (n — 2)-forms 
such that 

fsa - a = dXs . 

Note that since fs is assumed to be smooth with respect to 6, we can take Xs to be smooth as well. 

Let dJ^ = o-Q — /o(cro) for the proper boundary o"o and define Uq = aoHJ^^ so that the boundary 
of J^^ is of the form 

dJ^+ = a+ - fo{a+) + 770 . (39) 



According to (22 ) 



Ss 



a 



a + 



Using the adapted deformations ipl'^ to map J-^ to Si^ and Ut, respectively, then gives 



£6 



[{rsTa-m*a]+ / , [{rsr\s-m*xs] 



(40) 



(41) 



Differentiate this relation with respect to 8 and use the definition (28) to obtain 

17 = 



d_ 
85 



(5=0 J£s 



(42) 



where La is the Lie derivative (63). 



Cartan's formula (64) implies that 



L^a = i^da + d{i^a) = + d{ii^a) , 
LaAo = iAdAo + c?(iAAo) = iA(/oa) - ^A" + d{ii^\o) , 



where we used da = and the exactness of /q. Substitution of these results into (|42|) yields 



d_ 
d6 



<5=0 J£s 



iA^ + N 



where 



iAtt + 



«A(/oa) 



il^a + / iAAo • 
hat 



(43) 



Consequently, in order to prove the theorem we must finally show that = 0. Using (39) for dT 



then (43) becomes 



N 



[«A(/oa) - /o(^Aa)] + / iAa+ / «aAo 

Jm J dot 



'Vo 



(44) 



«aAo • 
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Note that this expression for N is independent of the choice of a. In order words, if a is any 



form such that da = ^ and f^a — a = dXs then we may begin again using the new forms in the 
computation of (40), obtaining finahy (44) with a ^ a and Aq Aq. In particular, if /o is exact 
with respect to a, then it is also exact with respect to 5 = (/q)'^"^ a, for all A; G N. Using this new 



form in ( 44 ) gives 



N 



*(A-/*A) 



A)(/o)'« + 



iA(/o)'Ao . 
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Finally, note that is also independent of the choice of A: since (41) is independent of the 



choice of adapted deformations, (42 ) is also independent of the choice of A — this equation and the 



subsequent ones are valid for any deformation A that comes from a pair of adapted deformations 
ipl'^- In particular, we will use the deformation (37), so that A ^ A, which is normal to JF and 
thus A|^o = A^+ = since da^ C ??o C JT. Thus if we define y = A - /o*A then y G TJ^ by (30) 
and 

N = [ ivUofoi. 

Since the (n — l)-form a is regular when dim(ylo), dim(i3o) <^ ^^r^; assumed in (HI), Lem. [l] 
implies that as k ^ oo, — > 0; however, N is independent of k, therefore = 0. □ 



7 Examples 



Here we give two simple examples demonstrating that (21) gives the expected results for the lobe 
volume. In the first case, the lobes are simply given by the cross product of the lobes for an area 
preserving map with a circle. In the second, the map is the time T map of a nonautonomous flow. 



7.1 Semidirect Product with a Twist Map 

An area-preserving, twist map {X, Y) = g(x, y) can be obtained from an implicit generating function 
G(x, X) through the equations 

y = d2G{x, X), y = -diG{x, X) . (45) 

provided that G satisfies the twist condition did2G ^ 0. A map generated in this way is also exact 
in the sense of Indeed we may choose u = ydx and the zero-form 

S{x,y) = G{x,X{x,y)) , (46) 

for in this case 

g*z/ — u = YdX — ydx , 
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and 

dS{x, y) = diG{x, X)dx + d2G{x, X)dX , 



which are equivalent to ( 45 ) . 

The map g can be extended to an exact volume-preserving map on M = X S by introducing 
an angle G S = M/ (27rZ), and some appropriate dynamics. A simple case is the semidirect product 

{X, y, 9) = /(x, y, 9) = {g{x, y),e + p{x, y)) , (47) 

where p is the local rotation number. Now / is an exact volume-preserving map with the volume 
form 0, = dx A dy A dO. For example, using the two form 

a = d{So /-I) AdO- ydx A dd , 

so that Q = da, and dS = YdX — ydx gives 

f*a-a = dS AdQ - YdX A dS - d {S o f-^) Ad9 + ydx A dO 
= -d{So /-I) Ade-y {dyp) dxAdy. 

Therefore, / obeys ^ with the one form 

\=-{Sor^)de-<i>{x,y)ydy. (48) 

where (/> is a function such that dx(t) = dyp. 

Because of the semidirect product structure of (47), if a is a saddle fixed point of g, then 
A = a X S is a normally hyperbolic invariant circle for /. Similarly if 7^ is a resonance zone for g 
(recall for example Fig. [l| then = 7^ x S is a resonance zone for /. Moreover, if £ is the exit set 
for IZ that is delineated by a pair of primary intersection points p and g, then £ = £ y.'B'is the exit 
set for the three-dimensional resonance zone of /. Moreover, since for each (x, y) the map Q = 6 + p 
is a homeomorphism of the circle S, the primary intersection manifold of the fundamental domains 
for / is the submanifold 

V = Vp-Vg (49) 
where rjp = {p} x S and rjq = {q} x S, and p, q are in the primary homoclinic orbits for g, recall 
Fig.[T] 

Thus we can apply Thm. |4] to compute the volume £ x S 

voi(o = e/ u'r^ = Y.\l ^- 1 

where A is the one form (48). These sums converge without the assumption that A is regular 
because f^{r]p^q) A as k ^ ±00. Noting that f^{rip) = g'^ip) x S, and that dy vanishes on the 
orbits of the primary intersections, the integrals above can be easily computed. For example. 



A = -Sig'-'ip)) I d9 = -2TrS{g'~\p)) 

fHvp) 
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Hence the volume becomes 

Vol(^) = 2^ {Sig\q)) - S{g\p))) . (50) 

which 2tt times the difference between the actions of the orbits of q and p under the symplectic 
map g. The action difference is exactly the lobe area for an area-preserving map as shown in 
|MMP87l lEasQlj , so that we have derived the obvious formula 

Yol{£) = 27r area(£:) . 

As an explicit example consider the generator 

G{x,X) = ^{X-hix)f , 

for a diffeomorphism /i : M ^ M. This generator satisfies the twist condition since did2G = 
—h'{x) 7^ by assumption. The generated map is, explicitly, given by 



= (X.r) = (/.W + 4^,^ I . (51) 



For this case, the generator (46) is 

^(-■■^'=Ks^)'=^'"- '''' 

Note that the line y = is invariant under g, and the dynamics on this line is simply x ^ h{x). 
Thus if h has exactly two fixed points, say xi < X2, then g has saddle fixed points at {xi,0). Two 
of the manifolds of these saddles simply lie on the j;-axis; for example, if < h'{xi) < 1 < h'{x2) 
then W^(xi,0) = (— 00,2:2) and VF^(x2,0) = (xi,cxd). The other two manifolds, M^"(xi,0) and 
VF'^(x2,0), may also intersect |Tab95] . and give rise to a nontrivial resonance zone, see for example 
Fig. [7j For this example, there are exactly two primary intersection orbits in the upper-half plane, 
labelled p and q in the figure. 



The three-dimensional map / defined through (47) has an invariant plane {y = 0} and, when 



h has two fixed points, has a pair of hyperbolic invariant circles 

Ci = {[xi,^,e):e e^} , i = l,2. 

The set S = {(x,O,0) : xi < x < X2} = W^{Ci) H W^{C2) is a heteroclinic manifold that has zero 
flux. The manifolds W^iCi) and W^{C2) can be used to form the second, leaky boundary of a 
resonance zone. Indeed, if p G W^{Ci) n W^{C2) is a primary intersection point for g^ recall Fig. [7] 
then the circle r/p = {p} x S is a proper boundary for both VF"(Ci) and W^{C2)- We can then form 
a resonance zone IZ with the boundary 



^^ = ^uT^;(^^)(Ci)uTy;(^^)(C: 



-2 
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Figure 7: Phase space of the map (51) for h{x 



^ cos(27rx) with a = 0.85. Here the saddle 



points lie at xi = —X2 = —0.25 and there is an elliptic fixed point at (0, ^). 



Notice that in this case, no "caps" are needed. If, as in Fig. [7j g has exactly one more primary 
intersection point, q, the set of primary heteroclinic intersections between Ci and C2 is given by 



(49). Finally, by (50) and (52), the volume of the exit lobe is 



fcez 

7.2 Nonautonomous Hamiltonian flow 

The flow of an incompressible vector field is volume- preserving. If the vector field is exact, in the 
sense we state below, then its flow will be exact as well. In this section we compute the volume of 
lobes that are obtained from the time-T map of such a flow. For this case, our results reduce to 
those of |Mac94j . 

Recall that an incompressible vector field X satisfies Lx^ = (V • X)Q = 0, where Lx is the 



Lie derivative. By (64), X is incompressible if the form ix^ is closed. Consequently, it is natural 



to say that X is exact-incompressible when ix^ is exact. Indeed this implies that the flow of the 
vector field is exact volume-preserving. 



Proposition 7 ( |LM08| ). Suppose X : M — > TM is a vector field with complete flow ipt and 
0, = da is an exact volume form on M . If ix^ is exact, then Lxa = d(3 for some n — 2 form (3, 
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and the flow is exact volume preserving, ip^a — a = dXt with the (n — 2)-form 

Xt= [ ip*Jdr, (53) 



for each t G M. 

As an example, consider the case of a 1 + ^ degree-of-freedom Hamiltonian flow, generated by 
a function H{x,y,t). On the extended phase space z = (x,y,0) £ H generates the vector 
field 

For the volume form Q = dx Ady AdO we can choose a = —ydx A dO so that da = 0,. Hence, 

ix^ = dH AdO + dx Ady , 

ixoi = —y——dO + ydx . 
oy 

These imply 

(dH\ 
H — w-;— ] AdO , 
dy J 

so that we can define P = —CdO where 

C = ydyH - H 

is the (phase space) Lagrangian. Assuming that the flow (ft of the vector field X is complete, then 
Prop. [7] implies that it is exact volume preserving with 

\t = -( f Coif^dT^de , (55) 



by (53). Thus Xt{X) is the (negative of the) action of the orbit segment from r = to t. 

Now suppose that H is 1— periodic in its last coordinate: H{x,y,6 + 1) = H{x,y,0), and let 
M = X T, where T = R/Z. 

The Poincare return map to the section 6 = is 

P{x,y) = ipi{x,y,0) . 

We will assume that P has a saddle fixed point, a, whose stable and unstable manifolds intersect 
and have exactly two transversal, primary intersection orbits, e.g., the orbits of p and q as shown 
in Fig. [!} This implies that the vector field (54) has a hyperbolic invariant circle A = {ipt{a,0) : 
< t < 1} whose two-dimensional stable and unstable manifolds have two, primary intersection 
orbits, (pt{q) and (pt{p) as sketched in Fig. [sl 



31 



WXA) 





Figure 8: Resonance zone for a periodically time dependent flow with a hyperbolic invariant circle 
A and two primary intersection orbits x^- 



By reducing to the Poincare section, the resonance zone and its exit and incoming sets for this 
three-dimensional system can be obtained purely by considering the area-preserving map P. In this 
case, we can compute the area of the exit lobe using the standard theory of |MMP84] . to obtain 

area(£:) = AW (56) 

where AW is the difference between the actions of the two primary homoclinic orbits: 

/oo 
{CiMQ))-^{Mp)))dT (57) 

Note that since P is the time-one map of X, area(£') is the area per unit time that escapes from 
the resonance zone. 

To make a test case for Thm. [4j we now consider the three-dimensional map 

f{x,y,e) = ipT{x,y,e) 

for T Z. By Prop. [7j / is an exact volume-preserving diffeomorphism of M with one form A-p 



given by (55 ). Under the above assumptions, j4 is a normally hyperbolic invariant circle of the map 



/, and / satisfies (H1)-(H3). 
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To construct a resonance zone, TZ, we must introduce caps: since the orbits of p and q are 
homotopic to A, they cross any proper loops a and 7. As usual we select a pair of proper loops 7 
and a, their corresponding fundamental domains U and S, and a cap C obeying (H4)-(H6). We 
let £ denote the three-dimensional exit lobe for this system. Using the form (55), since ^^d9 = dO 
we have 

(/')*At 



C o ip^dT dO 



»(fc+l)T 



C o {p^dT dO 



(58) 



The set 77 C V{A,A) used in Thm. |4] corresponds to the primary intersections on the fundamental 
domains, thus, = Vp ~ Vg where 

i]p = {^t{p) :t£R}nU = Ws{p*) : < s < T} , 
rjq = {yptiq) :teR}nU = Ws{q*) : < s < T} , 

and {p* , q*} = ^ r\ a are points on the orbits of p and q. 

Theorem [4] implies that, in order to compute the volume, we must integrate the one form given 
in (58) over r] and sum over /c G Z. Using the temporal parametrization for ijp gives d9\rjp 



ds, 



and 



(/')*At 



Vv 




{k+l)T 



C{ipt+s{p*))dt ds 



-T 



(k+l)T 



kT 



C{MP*))dr 



where we defined a new integration variable t = t + s and used the periodicity of C in time. The 



volume of the exit lobe for / is given by (21), and the sum becomes a single integral: 

/•CO 

Yol{£) = -T / {C{Mp*)) - ^Mq*))) dt . 



This integral is geometrically convergent since the orbits of p* and q* are bi-asymptotic to the 
hyperbolic circle A. Moreover, we can replace the initial conditions with p, and q, and thus obtain 



Vol(£:) = TAW , 



(59) 



where AT^ is given by (57). Note that Yol{£) is the volume that exits from the three-dimensional 
resonance per step of the map /, that is per T units of time. Thus, (59) is exactly what is expected 
from (56). 



8 Conclusion and future research 

The computation of lobe volumes is a first step toward developing a theory of transport for volume- 
preserving maps. These maps appear to have all sticky regions surrounding invariant tori and 
algebraic decay of exit and transit time distributions, complications that are familiar from the 
study of area-preserving maps |Mei92] . Since volume-preserving dynamics pertains to the motion 
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of Lagrangian tracers in incompressible fluids, a theory of transport should prove useful for the 
understanding laminar mixing and for designing an optimal mixer |Bal05j . 



To apply our result (21) to a general map we must compute the primary intersection V{A,B) 
between a pair of invariant sets A and B. Finding this set is easiest when the map / is reversible, 
R o f = o and when R{A) = B for the reversor R. In this case, it is easy to see that if 
the manifold M^'^(A) intersects the flxed set of the reversor, Fix(i?), the intersection point must be 
a point in the primary intersection set V. Thus this point can be used as a starting point for a 
continuation method to obtain rj. The intersection with Fix(/ o R) will give a point on a second 
component of r]. An suitable example to study is the quadratic map of |LM98] , which is reversible 
in special cases. We hope to report such computations in a future paper. 

A Some Notation 

Here we set out our notation, which follows e.g. |AM78j . If a is a fc-form and Vi,V2, ■ ■ - Vk are 
vector flelds, then the pullback, /*, of a diffeomorphism / is deflned by 

{raUVi,V2,...,Vk)=af^,j{Df{x)V,{x),...,Df{x)Vk{x)) . (60) 

The pullback can be applied to a vector fleld V as well: 

{rV){x) = {Df{x))-^V{f\x)) . (61) 

The pushforward operator is deflned as 

/* = {f-'T ■ 

The inner product of a with V is deflned as the {k — l)-form 

iya = a{V, •,...,•) . (62) 

Suppose that ipt is the {C^) flow of a vector fleld V, so that (po{x) = x, and d/dtipt{x) = V{(pt{x)). 
Then the Lie derivative with respect to V is the differential operator deflned by 



r - ^ 



ipt- (63) 

i=0 



where • is any tensor. The key identity for the derivative is Cartan's magic formula: 

Lv = iv{d) + d{iv) (64) 
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